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Abstract 

We find that the exponential operator V = exp [iA (Q1P2 + Q2P3 + • • • + Q„_i-P n + Q n fi)] , 
Qi, Pi are respectively the coordinate and momentum operators, is an n-mode squeezing operator 
which engenders standard squeezing. By virtue of the technique of integration within an ordered 
product of operators we derive V's normally ordered expansion and obtain the n-mode squeezed 
vacuum states, its Wigner function is calculated by using the Weyl ordering invariance under 
similar transformations. 

PACS 03.65.-w — Quantum mechanics 
PACS 42.50.-p — Quantum optics 



1 Introduction 

Quantum entanglement is a weird, remarkable feature of quantum mechanics though it implies 
intricacy. In recent years, various entangled states have brought considerable attention and interests 
of physicists because of their potential uses in quantum communication [TJ [5] . Among them the 
two-mode squeezed state exhibits quantum entanglement between the idle-mode and the signal- 
mode in a frequency domain manifestly, and is a typical entangled state of continuous variable. 
Theoretically, the two-mode squeezed state is constructed by the two-mode squeezing operator S = 
exp[A(aia2 — ala^)] [31 HI IS] acting on the two-mode vacuum state |00), 



S |00) = sechAexp — a 1 a 2 tanhA |00) , 



(1) 



where A is a squeezing parameter, the disentangling of S can be obtained by using SU(1,1) Lie 
algebra, [aia2, a\a\\ — a i a i + a \ a 2 + 1) or by using the entangled state representation \r) = r)i + iT] 2 ) 

mm 



1 77) = exp 
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100), 



(2) 



I -77) is the common eigenvector of two particles' relative position (Qi — Q2) and the tota momentum 
(Pi + P2), obeys the eigenvector equation, (Q\ - Q 2 ) \rf) = V^m Iv) > (-Pi + -P2) = \v) = V^m \rf) , 
and the orthonormal-complete relation 
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(3) 
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because the two- mode squeezing operator has its natural representation in (77 1 basis 
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b = cxp A \ a1a2 — a[a 2 J = J 



-),M = e A , (4) 
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The proof of Eq. Q is proceeded by virtue of the technique of integration within an ordered product 
(IWOP) of operators 0EECEU] 

7T/X 




(ai - a^j (a[ - a-A 



^ — a 1 a 2 tanh Ag(a 1 ai +a 2 Q2 + 1) In scchA^aia2 tanh A _ /g\ 



Eq. ([4]) confirms that the two-mode squeezed state itself is an entangled state which entangles the 
idle mode and signal mode as an outcome of a parametric-down conversion process The \rj) 

state was constructed in Ref. [5J [7] according to the idea of Einstein, Podolsky and Rosen in their 
argument that quantum mechanics is incomplete [12j . 

Using the relation between bosonic operators and the coordinate Qi, momentum Pj, Qi = (a, + 
a\)/\/2, Pi = (ai — a|)/(\/2i), and introducing the two-mode quadrature operators of light field as 
in Ref. [4j, x\ = (Qi + Q2V2, x-i — (Pi + P2V2, the variances of x\ and X2 in the state S |00) are 
in the standard form 

(00| S^xjS |00> = i e - 2A , (00| S^xjS |00> = ie 2A , (6) 

thus we get the standard squeezing for the two quadrature: On the other 

hand, the two-mode squeezing operator can also be recast into the form S — exp [iA (Q1P2 + Q2P1)] ■ 
Then an interesting question naturally rises: what is the property of the n-mode operator 

V = exp [iA (Q X P 2 + Q2P3 + ■■■ + Qn-lPn + QnPl)] , (7) 

and is it a squeezing operator which can engenders the standard squeezing for ?i-mode quadratures? 
What is the normally ordered expansion of V and what is the state V |0) (|0) is the n-mode vacuum 
state)? In this work we shall study V in detail. But how to disentangling the exponential of VI 
Since all terms of the set QiPj+i (i = 1 • • • n) do not make up a closed Lie algebra, the problem of 
what is V's the normally ordered form seems difficult. Thus we appeal to the IWOP technique to 
solve this problem. Our work is arranged in this way: firstly we use the IWOP technique to derive 
the normally ordered expansion of V and obtain the explicit form of V |0); then we examine the 
variances of the n-mode quadrature operators in the state V |0), we find that V just causes standard 
squeezing. Thus V is a squeezing operator. The Wigner function of V |0) is calculated by using the 
Weyl ordering invariance under similar transformations. Some examples are discussed in the last 
section. 



2 The normal product form of V 



In order to disentangle operator V, let A be 

/ 1 
1 




\1 
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o / 



(8) 



then V in Q is compactly expressed as 

V = exp 



(9) 



Using the Baker-Hausdorff formula, e^Se^ = B+ [A, B] + i [A, [A, B]] + ± [A, [A, [A, B]]] + ■ ■ ■ ,we 
have (here and henceforth the repeated indices represent the Einstein summation notation) 



V~ l Q k V = Q k -XQ i A ik + -iX 2 [Q i A ij P j ,QiA lk ]^ 

= Qi(e- XA ) ik = {e- xA ) kl Q u 
V- l P k V = P k + XA ki P i + ^iX 2 [A ki P j ,Q l A lm P m ] 



(10) 



(11) 



From Eq. fTUl) we see that when V acts on the n-mode coordinate eigenstate \q) , where q 
(qi, Qii • • • ! Qn)i it squeezes \q) in the way of 



V \g) = |A| 1/2 \Aq), A = e- XA , |A| = det A. 
Thus V has the representation on the coordinate (<f| basis 

V = J d n qV \q) (q\ = |A| 1/2 / d n q \Aq) (q\ , = V~\ 

since J d n q \q) (q\ = 1. Using the expression of eigenstate \q) in Fock space 

\q) = 7T-™/ 4 : exp[-±<fg + - §aV] |0) , 

a f = (a\,al, ■ ■ ■ ,o+), 
and |0) (0| =: exp[— &tat] : , we can put V into the normal ordering form , 

V = vr-™/ 2 \A\ 1/2 J d n q: exp[-^f(l + AA)g + VTq(ka) + a) 

— — (da + d^a^) — a^a] : . 



(12) 



(13) 



(14) 



(15) 



To compute the integration in Eq. lTSl) by virtue of the IWOP technique, we use the mathematical 
formula 

"1. 



J cTxexp[-xFx + xv] = vr™ /2 (det F)~ 1/2 exp 



-vF^v 



(16) 



then we derive 



V 



det A 



1/2 



exp 



-a* (AAT^A- /V 



det TV 

: exp [a f (AA^ 1 - /) a] : exp 



(17) 



where AT = (1 + AA)/2. Eq. (fT"7|) is just the normal product form of V. 



3 The squeezing property of V |0) 

Operating V on the n-mode vacuum state |0) , we obtain the squeezed vacuum state 



V\Q) 



det A 
detiV 



1/2 



exp 



-at [an- 1 



A - I ) a f 



(18) 



Now we evaluate the variances of the n-mode quadratures. The quadratures in the n-mode case are 
defined as 



A', 



v 2n . 

2 /v2\ /xr\2 



(19) 



obeying [X\, X%\ = |. Their variances are (AXj) = (Xf} — (X~j) , i = 1,2. Noting the expectation 
values of X\ and X2 in the state V J 0) , (Xi) = (X2) = 0, and using Eqs. (fTU|) and (|TT|) we see that 
the variances are 



. n n 

(ax,) 2 = (o\v- 1 xfv\o) = -(o\v- 1 Y / Q*Y,Q3 v \° 
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E( e ~ AA )«( e ~ Aj[ )^ - 4^E( Xa )«> 
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similarly we have 



1 " 

(AI2) 2 = <0| F-^fy |0> = -^(AA)r.i, 



(20) 



(21) 



Eqs. ([20]) - ([2T|) are the quadrature variance formula in the transformed vacuum state acted by the 

n 

operator exp[iA Qi-^ijPj]- By observing that A in (|8|) is a cyclic matrix, we see 



n 



2'n, 



* j 



then using = AA, so AA = e A ( A +- 4 ) i a symmetric matrix, we have 



1=0 (j 



= n >^ v ,, - 2' = ne 



(-A)< 



)( _ „„-2A 



and 



it then follows 



E (^A)rj = 



1 ~ 

( A ^) 2 - t-E( aa ) 



,2A 



-2,\ 



2A 



(AX 2 ) 2 = ^E(AA),^V- 



(22) 



(23) 
(24) 

(25) 
(26) 



This leads to AX"i • AX2 = j, which shows that V is a correct n-mode squeezing operator for the 
n-mode quadratures in Eg. (|19[) and produces the standard squeezing similar to Eq. ([6]). 



4 The Wigner function of V |0) 



Wigner distribution functions |T3l [141 E] 01 quantum states are widely studied in quantum statistics 
and quantum optics. Now we derive the expression of the Wigner function of V |0) . Here we take a 
new method to do it. Recalling that in Ref . [THl El HHj we have introduced the Weyl ordering form 
of single- mode Wigner operator A (q,p), 



A, ( qi , Pl )= - J(«i-Oi)*(pi-i\) ', 



(27) 



its normal ordering form is 



= ~: exp - {qi ~ Qif - (pi - Pif 

7T L 



(28) 



where the symbols : : and ; ; denote the normal ordering and the Weyl ordering, respectively. Note 
that the order of Bose operators a\ and a\ within a normally ordered product and a Weyl ordered 
product can be permuted. That is to say, even though [ai,aj] = 1, we can have : a\a\: =: a\a\ : 



and ; a\a\ ; = ; a[a\ \ . The Weyl ordering has a remarkable property, i.e., the order-invariance of Weyl 
ordered operators under similar transformations [161 117[ 118) , which means 



U' (ooo) 'J/" 1 = '[/(ooo)?/- 1 ' 



(29) 



as if the "fence" ; ;did not exist. 

For n-mode case, the Weyl ordering form of the Wigner operator is 

A„ (q,f)= : d(q-Q)s(p-P 



(30) 



where Q — (Qi, Q2, • ■ • , Q n ) an d P — (Pi> P21 • • • ,Pn)- Then according to the Weyl ordering 
invariance under similar transformations and Eqs. (|10[) and (jlip we have 



V~ 1 A n (q,p)V 



V- i: s(q-Q^s(p-P) \ V 
8 (q k - (e- xA ) M Q t ) S ( Pk - (e rA ) kl P t ) [ 



6 e rA q 



= A 



Q)s(c 



-rA 



p- P 



Thus using Eqs.pT]) and d3TJ) the Wigner function of V |0) is 

(0|^- 1 A„ (q,f)V\0) 
= \ (0| : eM~(e rA q- Qf - {^ rA P~ P?\ ■ |0 



eM~(e rA q) 2 



(e- rA p) 2 } 



exp 



~ ex P 

7T™ 



-qe rA e rA q - pe~ rA e~ rA p 



-q(AA) q-pAAp 



(31) 



(32) 



From Eq.((32|l we see that once the explicit expression of AA = exp[— X(A + A)] is deduced, the 
Wigner function of V |0) can be calculated. 



5 Some examples of calculating the Wigner function 



Taking n — 2 as an example, V n=2 is the usual two- mode squeezing operator. The matrix A 
j q ) . it t hen Follows t hat 



and 



AA = e -M*M) 



AA 



cosh2A — sinh2A 
— sinh 2A cosh 2A 



cosh 2A sinh 2A 
sinh2A cosh2A 



(33) 
(34) 



Substituting Eqs.fO and flM|) into Eq.(|3"2). we have 



(00\V- 1 A 2 (q,p)V \00) = ^exp |-2 (a*a% + a x a 2 ) sinh2A - 2 (joi| 2 + \a 2 \ 2 j cosh2Aj , (35) 

where on — (qi + ipi) , (i = 1,2.). Eq. (|35|) is just the Wigner function of the usual two-mode 
squeezing vacuum state. For n — 3, we have 



u v v i 2 1 I 

AA= | v u v , u = -e x + -e- 2X , « = - ( e ~ 2A 

V V u '* *' 



and 



^AA^ is obtained by replacing A with —A in A A. By using Eq. 

3 



(36) 

the Wigner function is 



(0\V" 1 A 3 (q,p)V\0) 



■ exp 



- (cosh 2 A + 2 cosh A) |a< 
exp < — - (sinh 2 A — 2 sinh A) a 2 

- [(cosh 2A — cosh A) a^a* + (sinh A + sinh 2A) ajjoy] + c.c (<87) 



i>i=i 



For n — 4 case we have (see the Appendix) 



AA 



U W V w 



y w' v' w' u' j 



(38) 



where vl — cosh 2 A, v' — sinh 2 A, w' = — sinh A cosh A. Then substituting Eq. ([55|) into Eq. (j52"|) we 
obtain 



1 



(0|y _1 A 4 (q,p)V\0) = — exp<^ -2 cosh 2 A 



E 



(M + M*) tanh 2 A + (R* + R) tanh A 



(39) 

where M = a±a 3 + a 2 a* Al R — oi\u 2 + u\u^ + a 2 a 3 + 0:30/4. This form differs evidently from the 
Wigner function of the direct-product of usual two two-mode squeezed states' Wigner functions (I35[) . 
In addition, using Eq. (f38|) we can check Eqs. (f25|) and (|26| . Further, using Eq. ([38|) we have 



I 2 tanh A tanh A A 

tanh A 2 tanh A 

tanh A 2 tanh A 

\ tanh A tanh A 2 J 



det N = cosh A. 



(40) 



Then substituting Eqs. (|4"0)) and (A. 4) into Eq.JTH]) yields the four-mode squeezed state, 



V |0000) = sechAexp 



4) (4 + 4) 



tanh A 



|0000) , 



(41) 



from which one can see that the four-mode squeezed state is not the same as the direct product of 
two two-mode squeezed states in Eq.(TT]). 

In sum, by virtue of the IWOP technique, we have introduced a kind of an n-mode squeezing 
operator V = exp [iA (Q1P2 + Q2P3 + • • ■ + Qn-i-fn + QnPi)], which engenders standard squeezing 
for the n-mode quadratures. We have derived V's normally ordered expansion and obtained the 
expression of n-mode squeezed vacuum states and evaluated its Wigner function with the aid of the 
Weyl ordering invariance under similar transformations. 

Appendix: Derivation of Eq.(|38p 

For the completeness of this paper, here we derive analytically Eq. (j3"8)) . Noticing, for the case 
of n = 4, A 4 = I, I is the 4x4 unit matrix, from the Cayley-Hamilton theorem we know that the 
expanding form of exp(— rA) must be 



A = cxp(-A^) = c (A)J + a(\)A + c 2 (X)A 2 + c 3 (X)A 3 . 
To determine Cj(X) , we take A being e 1 ^/ 2 )" (j = 0, 1, 2, 3) respectively, then we have 
exp(-A) = c (A) + ci(A) + c 2 (A) + c 3 (A), 

exp(-Ae i ( 1 / 2 )") = cq(A) + ci(A)e i ( 1 / 2 ^ + c 2 (A)e iT + c 3 {X)e^ 3 ' 2 ^ , 
exp(-Ae i7T ) = co (A) + ci(A)e i7r + c 2 (A)e i27r + c 3 (A)e i37r , 
exp(-Xe^ 3 ^) = c (A) + c 1 (A)e i ( 3 / 2 )" + c 2 {X)e^/ 2 ^ + c 3 (X)e i( - 9 / 2 ^ . 

Its solution is 



(A.l) 



(A.2) 



It follows that 



and 



A 



co (A) = h (cosh A + cos A) 
ci (A) = J (— sinh A — sin A) 
c 2 (A) = \ (cosh A — cos A) 
C 3(A) = I (— sinh A + sin A) 



/ c c 3 c 2 Ci \ 

ci c c 3 c 2 

c 2 Cl Co c 3 

y C 3 C 2 Cl Co J 



, det A = 1, 



(A.3) 



(A.4) 



AA = [c (A)7 + Cl (X)A + c 2 (X)A 2 + c 3 (A)A 3 ] • c (A)/ + a{X)A + c 2 (A)i 2 + c 3 (A)i 



/ 2 cosh 2 A -sinh2A 



sinh 2A 2 cosh A 



2 sinli A 
— sinh 2A 



2 sinli A 
— sinh 2 A 



sinh 2 A 2 cosh A 



-sinh2A \ 
2 sinh 2 A 
— sinh 2 A 



(A.5) 



2 sinh 2 A -sinh2A 2 cosh 2 A / 



this is just Eq. (|38|) . 
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